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ABSTRACT 

Using U'- through i^ s -band imaging data in the GOODS-South field, we construct a large, complete 
sample of 275 "extremely red objects" (EROs; K s < 22.0, R-K s > 3.35; AB), all with deep HST/ ACS 
imaging in -B435, P606, *775j and zsso, and well-calibrated photometric redshifts. Quantitative concen- 
tration and asymmetry measurements fail to separate EROs into distinct morphological classes. We 
therefore visually classify the morphologies of all EROs into four broad types: "Early" (elliptical- like) , 
"Late" (disk galaxies), "Irregular" and "Other" (chain galaxies and low surface brightness galaxies), 
and calculate their relative fractions and comoving space densities. For a broad range of limiting 
magnitudes and color thresholds, the relative number of early- type EROs is approximately constant 
at 33-44%, and the comoving space densities of Early- and Late- type EROs are comparable. 
Mean rest-frame spectral energy distributions (SEDs) at A rcst « 0.1-1.2/im are constructed for all 
EROs. The SEDs are extremely similar in their range of shapes, independent of morphological type. 
The implication is that any differences between the broad-band SEDs of Early-type EROs and the 
other types are relatively subtle, and there is no robust way of photometrically distinguishing between 
different morphological types with usual optical/near- infrared photometry. 

Subject headings: surveys - galaxies: evolution - galaxies: fundamental parameters - galaxies: high- 
redshift - galaxies: stellar content - infrared: galaxies 



1. introduction 

The newest challenges to galaxy formation models lie 
in predicting or explaining the state of galaxy assembly 
when the universe was in its middle ages, at only several 
billion years old. In field populations at these redshifts, 
z ~ 1-2, it is necessary to sample a large enough vol- 
ume at sufficient depth to overcome biases due to cosmic 
variance. This is particularly true for strongly clustered 
galaxies, such as the population of "extremely red ob- 
jects" (EROs). These have comoving correlation scales 
of ~ lOfe^nMpc at z > 1 l IDaddi 1 et a.U TMM HI 
iFirth et al J 120021: iRoche et all 12002(1 . and therefore re- 
quire surveys of solid angle at least lOOarcmin 2 to have 
a variance less than ~ 20% l)Somerville et al.ll2003ft . 

The spectral energy distributions (SEDs) of galaxies 

1 Based on observations taken with the NASA/ESA Hubble 
Space Telescope, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc. (AURA) under NASA contract 
NAS 5-26555. 

2 Based on observations collected at the European Southern 
Observatory, Chile (ESO Programmes 164.O-0561, 168.A-0485, 
169.A-0725, 170.A-0788, 267.A-5729) 

3 Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218; leonidas, casertano, med, ferguson, 
mauro, koekemoe, somerville@stsci.edu 

4 California Institute of Technology, Mail Stop 105-24, Pasadena, 
CA 91109; cc@astro.caltech.edu 

5 Jet Propulsion Laboratory, California Institute of Technology, 
Mail Stop 169-327 (PE) & 169-506 (DS), Pasadena, CA 91109; 
prmeOkromos . jpl .nasa. gov , sternOzwolf kinder . jpl . nasa .gov 

6 Department of Physics and Astronomy, The Johns Hopkins 
University, 3400 North Charles Street, Baltimore, MD 21218; 
nagroginOstsci . edu 

' Steward Observatory, The University of Arizona, 933 North 
Cherry Avenue, Tucson, AZ 85721: papovichOas.arizona.edu 

8 European Southern Observatory, Karl-Schwarzschild-Strasse 2, 
D-85748 Garching bei Miinchen, Germany; renzini0eso.org 



are products of the amalgam of stellar mass functions, 
star formation histories, dust, and geometry, often with 
degeneracies (such as in age and metallicity) compli- 
cating interpretation. EROs, by their simple optical- 
minus-infrared color selection, embody this degeneracy 
extremely well, and are known to consist of at least two 
distinct classes of galaxies: old, and dusty/star-forming. 
Surveys in other wavelengths have post facto uncov- 
ered EROs, including r adio iWillott et~all I2001D . sub- 
mm l|Smail et all 120021: iWehner et all 12002ft . and both 
hard- and soft-X -ray surveys IjAlexander etUl 12001 
IBrusa et al.ll200 2L Each class (and subclass) of EROs 
offers a potentially important datum for further unfold- 
ing and challenging galaxy formation models. Currently, 
predictions of CDM-based semi-analytic models do not 
compare w ell against the numbers and statistics of EROs 
as a whole ijCimatti et al.l2002rilFirth et alJ2002]) . Hav- 
ing the relative fractions and space densities of, mini- 
mally, morphologic ally-determined early-typ e EROs ver- 
sus all EROs (e.g. lYan fc ThonTDSonll2003|L should im- 
prove the situation for comparisons against models. 

Currently, there is a popular working hypothesis that 
some (large) fracti on of EROs (those similar to the 
z = 1.552 53W091. ISpinrad et al.lll997|) ar e the ances- 
tors of present-day elliptical galaxies (e.g. IDaddi et alJ 
I2000ri IMoustakasfc Somervilld 12001 [MS02]; and ref- 
erences therein). Conversely, many EROs demonstra- 
bly host high ly obscured s tarbursts (e.g. the z = 
1.44 HR10: iGraham fc Devi Il99^ which could be an 
important component of the global lum inosity density 
l|Cimatti et al.ll2002al: lElbaz et al] 12002ft . A photomet- 
ric tec hnique has been proposed bv lPozzetti fc Mannuccil 
(2000) to make broad distinctions between passive and 
star-forming, dusty EROs. This is as yet not well-tested. 
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A systematic morphological and photometric study of a 
large, homogeneously-selected sample of EROs is there- 
fore timely. 

In this Letter we present the selection, morphologi- 
cal mix, and mean rest-frame SEDs of K s < 22.0 (AB) 
EROs , selected from the G OODS-South HST / ACS mo- 
saic i|Giavalisco et al.ll2003j) over 163arcmin . The sam- 
ple selection and their rcdshifts are described in §0 
followed by the morphological study in §0 In §0] we 
present the relative fractions, space densities, and av- 
erage rest-frame SEDs by morphological type. Conclu- 
sions are drawn in §0 We give magnitudes in the AB 
system 9 , and, where necessary, we use a flat cosmology 
with fi m = 1 — tt\ = 0.3 and Hq — 70 h?o kms -1 Mpc -1 . 
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Fig. 1. — The ERO rcdshift distributions N(z) for the full sample 
(R — K a > 3.35, K s < 22.0) and for each morphological type, as 
given in legend. The redshift distributions for Early- and Late- 
types are very similar, with z mc( j ~ 1.2, tailing to z ~ 1.5. The 
Irregular- and Other- types have a very broad N(z) distribution 
with 2 mct j PS 1.5, and dominate the sample at z ^ 1.5. Of the 
275 objects plotted here, 66 redshifts are spectroscopic and the 
remainder are photometric (see text). 



2. SAMPLE SELECTION AND REDSHIFTS 

The GOODS data are described in iGiavalisco et alJ 
(2003). Briefly, for the South ern field (overlapping with 
CDFS; IGiacconi et all 12002). an ACS mosaic covering 
163arcmin 2 in -B435, V§Q6, *775> and z^5o are supported by 
a broad set of ground-based optical and infrared imaging, 
U'UBVRIJHKs, from the VLT, NTT, and 2.2 m tele- 
scopes. The infrared data are from the SOFI and ISAAC 
instruments. For this work, we employ a PSF-matched 
SOFI ^-selected catalog from which the stars have been 
removed 10 . The SOFI data reach ~ 50% completeness 
at K s 2 2.8, and are still nearly entirely compl ete at 
K s « 22.0 ijGiavaiisco et al.ll2003UMov et al.ll2003|) . 

We consider several definitions of EROs, all of which 
are essentially complete for these survey depths: three 
color cuts (R — K s > 3.35, 3.65, and 4.35) for each of 
two magnitude limits (K s < 21.0 and 22.0). The largest 
("full") sample (i? — K s > 3.35, K s < 22.0) consists 
of 275 objects. The smallest sample {R — K s > 4.35, 
K s < 21.0) contains 18 objects. 

9 The most relevant conversions to Vega magnitudes and colors 
for this paper are: -Rvega = R — 0.206 and X 3 vega = K s — 1.859. 
Thus, R — K B > 3.35 in AB corresponds to (R — A" s )v C ga > 5.0. 

10 Unresolved sources are identified by a plot of the zg$o 
0.0625 arcscc-diameter surface-brightness vs isophotal magnitude, 
which provides a robust separation to zg$Q 26.5. 



Phot ometric re dshifts are computed using the method 
of lBenited l|2000D . These are well-calibrated to z SD0C w 
1.5, with cr[Az/(l + z spec )] = 0.11 (to R < 25.5; 
iMobasher etall l2003) using a total of 434 spectroscopic 
rcdshifts, including 202 from a VLT/FORS2 campaign 
over 2002-2003 (to be presented elsew here), and from 
the K20 survey ijCimatti et alJl2002bld) . This sample in- 
cludes 66 EROs with R - K s > 3.35 and K s < 22.0, 
for which cr[Az/(l + z S pcc)] = 0.05. In the analysis that 
follows, photometric redshifts have been replaced with 
spectroscopic redshifts where possible. There is one faint 
compact ERO with a formal z p hot = 6.4, which is a candi- 
date L-dwarf; for completeness we keep this object in the 
full sample. The ERO redshift distribution for the full 
sample is shown in Fig. ^ The N(z) has a fairly sharp 
edge at z « 0.95, a median value of z mo d ~ 1-2. The 
N(z) for the redder subsets are similar, but shifted to a 



higher median redshift, 



1.5 (for R-K s > 4.35). 
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Fig. 2. — Montage of 3-color Biz stamps, illustrating exam- 
ples of each morphological type, as described in the text. Two 
objects are given for each type (by column). The caption of each 
stamp gives the ID ("GS"=GOODS-South); the redshift (spectro- 
scopic for all but for GS 106 and GS3041, for which photometric 
redshifts are given); and a letter code reflecting observed spec- 
troscopic features — A=absorption; E=emission; C=both emission 
and absorption features. 



3. MORPHOLOGICAL CLASSIFICATION OF EROS 

The study of EROs has been largely motivated by 
their tantalizing connection with both the progenitors of 
present-day massive elliptical galaxies, and possibly, ul- 
traluminous infrared galaxies. Using the high-resolution 
ACS data, it should be possible to distinguish between 
these extremes morphologically. With the ACS z^so 
data, restframe morphologies longward of 4000A can 
be determined for redshifts z ^ 1.3, r educing possible 
morp hological fc-correction issues (e.g. IGiavalisco et all 
1996) for most EROs in this sample. The GOODS 
data reach a zssn l-c is ophotal limit of 27.3magarcsec -2 
ijGiavaiisco et all 12003). As the average magnitude of 
galaxies in the full sample is zsso ~ 23.0, for most cases 
the distinction between early- and non-early-types should 
be clear. At redshifts beyond z « 1.3, the 4000A-break 
begins to shift past the Z850 band, and morphological k- 
corrections may become important. This can be solved 
by also examining longer-wavelength imaging (e.g. the 
~ 0.4 arcsec-resolution ISAAC data in the K s -band); this 
will be done in future work. 

We began by exploring the galaxy morphologies 
with the concentration and asymmetry (CAS) param- 



Morphologies and SEDs of EROs in GOODS-South 



3 



eters of iConselicd iEM see also lAbraham et all Il996t 
iBershadv et alJl200nl) . We found that the parameter val- 
ues were fairly continuous and did not cleanly distinguish 
"classes," which reflects the morphological richness of all 
types of EROs. We therefore resorted to a visual clas- 
sification approach, and report the range and spread of 
the resulting by-type CAS values below. 

We adopt four broad morphological types: Early, 
Late, Irregular, and Other. The classification is based 
largely on the zsso data, with cross-examination of the 
-B435J ^606, and 2775 imaging for the colors and multi- 
wavelength structure of features. Since a primary goal 
of this work is to make a census of EROs that may be 
present-day elliptical progenitors, we use a conservative 
definition for Early types, excluding galaxies with any 
definite evidence of a prominent disk or dust-band. An 
isolated environment was not a requirement, so included 
are galaxies that are clearly interacting or merging. The 
Late-type category consists of galaxies with some reg- 
ular, disky component, independent of the bulge frac- 
tion - which can be quite large. Irregular-types are the 
Magellanic- type galaxies, whereas the Other type pre- 
dominantly includes low surface brightness galaxies and 
chain galaxies with multiple star-formation sites. Two 
examples of each type are given in Fig. [2J 

The classifications were done independently by four au- 
thors (SC, PE, RL, LAM), with the consolidation done 
by LAM. The overall agreement was very good, the great- 
est differences arising in the use of galaxies' central con- 
centration as a criterion for several cases of morphologi- 
cally Early-type galaxies. Since it became apparent that 
many EROs appear to be in dynamical disequilibrium, 
potentially softening the sharpness of cores, concentra- 
tion as a criterion was relaxed. 

For both Early- and Late- type EROs, the mean val- 
ues and typical dispersions of the CAS parameters, are 
C E « Cl « 3.1 ± 0.5, and A E w A h w 0.25 ± 0.1. The 
Irregular- and Other-type EROs are both less concen- 
trated (at C\ « Co ~ 2.5 ± 0.3) and more asymmetric 
(at Ai w A w 0.3 ±0.2). 

4. RESULTS 

The full ERO dataset comprises ~ 10% of all K s - 
selected galaxies at K s = 20-22. The relative fractions of 
EROs by morphological type are presented in Tablcnfor 
several limiting magnitudes and color limits. We find the 
relative fraction of Early-types does not depend strongly 
on the selection criteria, staying at 33 — 44% of the to- 
tal numbers. Late-types, however, dominate the popu- 
lation at brighter magnitudes (with 55% at K s < 21.0, 
R — K s > 3.35), and comprise much smaller fractions 
(23%) at the faintest, reddes t limit (see also ipsi 4iy-b and 
WFPC2-based ana lyses by lYan & Thompson! 120031 and 
iCilbank et alJ20T)3l which find similar relative fractions). 
The redshift distributions of Early- and Late-types are 
very similar (Fig.^). The combined fraction of Irregular- 
and Other-types increases dramatically with both limit- 
ing magnitude and color threshold; these galaxies also 
tend to have higher photometric redshifts, with a me- 
dian of z mo d ~ 1.5. Based on the observed redshift dis- 
tributions, comoving space densities are calculated and 
presented in Table [21 The volume probed by the full 
sample is ~ 4.4 x 10 5 Mpc 3 , for a total comoving 
space density of n = (6.40 ± 0.39) x 10" 4 h% Mpc" 3 , and 



the cosmic variance un certainty is therefore a cv « 15% 
llSomerville et a1jl20f)l. 

To explore whether there is a useful correspondence 
between morphological type and colors, we calculate the 
rest- frame SED of each ERO, normalized to f v = 1 (ar- 
bitrary units) at 5200 A. In Fig. we plot all data points 
for two groupings of morphological types, with the cor- 
responding mean "blue" and "red" SEDs, as described 
in the caption. The strong break around 4000 A fol- 
lows from the red-color selection of EROs. It is inter- 
esting that the range of SED shapes is extremely sim- 
ilar for all EROs. This is reflected also in the formal 
SED-fitting don e for ca lculating photometric redshifts in 
iMobasher et al.l l)2003|) . Even the reddest normal star- 
burst SEDs (drawn from lKinnev et al.ll996ft are too blue 
to fit EROs, so the EROs tend strongly to be best-fit 
by the earliest (elliptical) SED templates. As seen in 
the inset of Fig. |3J the empirical SED of the ultralu- 
minous infrared galaxy (ULIRG) Arp 220 (as given in 
lElbaz et al.1200 2'1 is sufficiently red to match ERO colors, 
and is coincidentally remarkably similar to the SED of an 
old, passively-evolved, unreddened system. Many of the 
Late-type galaxies have large bulges, and so their SEDs 
may naturally be dominated by old stellar populations. 
Even the Irregular and Other morphological types alone, 
which may arguably be candidates for ULIRGs, span the 
same range of SED shapes shown in Fig. U3 We conclude 
that in this wavelength-range, color information is insuf- 
ficient to generally distinguish between clearly distinct 
morphological types of EROs. 

5. CONCLUSIONS 

Samples of color-selected EROs consist of several dis- 
tinct galaxy types for all magnitude ranges and color- 
cuts. All morphological types are found over very 
similar redshift ranges, z mo d w 1.2, with the more 
morphologically-disturbed (Irregular and Other) types 
tending to somewhat higher redshifts. This may be due 
to a combination of a morphological /c-correction selec- 
tion effects (where at redshifts above z 1.3 galaxy 
bulges may be less apparent in the Z850 data); and a 
hint of a population of real, "train-wreck" galaxies, seen 
at redshifts closer to z « 2. With the present data, 
it is not possible to distinguish between these possibili- 
ties from the photometry alone, as seen by the grossly- 
similar rest- frame SEDs of all morphological types. This 
does suggest that those EROs that are perhaps "dusty 
starbursts," have SEDs that are more similar to local 
ULIRGs (such as Arp 220), rather than (mere) starburst 
galaxies. 

The comoving space density of morphologically Early- 
type EROs (tie ~ 2.5 x 10~ 4 /i 3 Mpc" 3 ) is quite compa- 
rable to present-day giant Ellipticals, n g E,z=o ~ 2.1 x 
10" 4 ^ Mpc" 3 (c.f. MS02). This suggests that the 
z w 1.2 Early- type galaxies identified within the ERO 
population may in fact correspond to the majority of the 
most massive ellipticals today, pointing to an epoch of 
assembly, as well as star-formation, at redshifts well be- 
yond z > 2. 

We acknowledge discussion at many stages of this work 
with D. Alexander and F. Bauer, and comments by L. 
Yan. The work of DS was carried out at JPL, Caltech, 
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Fig. 3. — The rest-frame SEDs for K s < 22, R — K s > 3.35 EROs, Early at left; the combined set of Late, Irregular, and Other types 
combined, at right. These are constructed from the broadband magnitudes, shifted to the restframe via photometric and spectroscopic 
redshifts. All SEDs are nor malized to f v = 1 at 5200 A, with the precise normalization value for each object determined via interpolation 
on the best-fitting SED from Mobash er et alj 12003). All points arc shown with their photometric uncertainty — some of the error bars are 
smaller than the points. At each wavelength bin, a bi-weight mean and dispersion are calculated; the blue and red "envelope" SEDs reflect 
approximately the mean bluest and reddest SEDs in eac h sample. These c an then be compared against two representative ERO templates, 
shown in the inset plot: the Arp 220 SED (as given in Elbaz ct al. 2002); and a 4.0 Gyr old Z = Zq simple stellar population from the 
v.2000 lBruzual fc Chariot 119931) models. 
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Table 1 

ERO FRACTIONS (%) BY MORPHOLOGICAL TYPE* 

K s < 21.0 K s < 22.0 

R — K 3 > 3.35 3.65 4~35~ 3.35 3.65 4.35 

~36 41 44 37 38 33~ 

55 52 39 38 32 23 

5.1 5.5 17 11 12 23 

2.5 1.4 0.0 15 17 20 

118 73 18 275 185 60 

*: The relative fractions are given in percent (%). 
* : Total number of EROs with shown criteria. 



early 
late 
irregular 
other 
total number^ 



Table 2 

ERO SPACE DENSITIES [xlO -4 hj Q MPC -3 ] BY MORPHOLOGICAL TYPE* 







K s < 21.0 






K s < 22.0 




R, — K B > 3.35 


R - K s > 3.65 


R — K B > 4.35 


R — K s > 3.3! 


5 R — K s > 3.65 


R, — K s > 4.35 


early 


1.10 ± 0.17 


0.67 ± 0.12 


0.09 ± 0.03 


2.49 ± 0.25 


1.59 ± 0.19 


0.29 ± 0.07 


late 


1.45 ± 0.18 


0.86 ± 0.14 


0.14 ± 0.05 


2.36 ± 0.23 


1.37 ± 0.18 


0.25 ± 0.07 


l+i+o f 


1.67 ± 0.19 


0.98 ± 0.15 


0.17 ± 0.05 


3.84 ± 0.29 


2.53 ± 0.24 


0.72 ± 0.11 


all 


2.94 ± 0.27 


1.64 ± 0.19 


0.26 ± 0.06 


6.40 ± 0.39 


4.11 ± 0.30 


1.03 ± 0.13 



*: The space density of EROs with specified R — K s and K s thresholds, in units of 10 4 h^ Mpc 3 , with Poisson uncertainties. 
Results for the combined numbers of Late-, Irregular-, and Other-type EROs. 



